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^ Abstract 

i } 

By  means  of  a systematic  approximation  based  on  the 
perturbation  theory,  from  t^e  three-dimensional,  non- 
linear elasticity  theory,  rive  partial  differential 
equations  are  obtained  for  the  non-linear  oscillations  of 
visco-elastic  olates.  The  Plate  can  undergo  a finite  trans- 
verse deflection.  The  rotatory  inertia,  shear  deformation  and 
effect  of  body  force  and  moment  components  are  taken  into 
account.  The  present  theory  is  the  generalization  of  Foppl- 
Karuan-Timoshenko  theory  of  large  static  deflection  on  the 
one  hand  and  Uf iyand-Mindlin  theory  of  small  oscillation 
theory  on  the  other  hand. 

Boundary  and  initial  conditions  are  discussed  and  an 
expression  for  the  strain  energy  is  given- 

1 . Introduction 

Foppi  obtained  two  simultaneous  non-linear  partial 
differential  equations  to  study  the  large  static  deflection 
of  plates  [1].  Von  Karman  and  3,  Timoshenko  reformulated  and 
studied  these  equations  for  various  plate  problems.  Many 
other  authors  used  these  equations,  now  known  as  FSopl-Karman- 
Timoshenko  theory.  Assumptions  or  this  theory  are  not 
systematic  and  thus  it  is  desiracie  to  know  whether  there  are 
other  terms  missing  from  these  equations  which  are  of  the 
same  order  of  magnitude  with  t.he  present  ones.  Also  for 
oynomic  problems  the  inertia  terms  must  oe  taken  ?iuo  account. 

/ 

The  effect  of  shear  c.ef ormation  m static  plate  problems 
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was  first  considered  by  £.  Reissner  [2].  Uflyand  [3]  and 
Mindlin  [4]  have  reformulated  small  oscillation  problems  of 
plates  and  took  into  account  the  effect  of  rotatory  inertia 
and  shear  deformation. 

Recently  for  the  general  formulation  of  shell  problems, 
many  papers  have  appeared,  most  of  which  are  in  tensor 
notation  or  written  in  intrinsic  coordinates  which  do  not 
lend  themselves  easily  to  application  [5,  6,  7]. 

The  highly  rigorous  nature  of  these  papers  leads  to  very 
high  complexity,  thus  making  the  solution  of  even  very  simple 
problems  prohibitive  or  impossible.  The  use  of  the  highly 
productive  method  of  engineers,  the  intuition,  is  thereby 
hindered. 

In  impact  problems,  where  the  transient  stage  is  important, 
a further  complication  must  be  included  into  the  formulation, 
namely  the  internal  damping  effect. 

With  these  views  in  mind,  the  reformulation  of  the  plate 
theory  is  made.  Consequently,  the  present  theory  includes  the 
effects  of:  (a)  finite  transverse  deflection,  (b)  shear 
deformation,  (c)  rotatory  inertia,  (d)  internal  damping, 

(e)  body  force  and  moments.  The  theory  is  based  on  the  following 
assumptions: 

1.  The  deformation  components  are  expandable  into  power 
series  of  non-dimensional  transverse  coordinates  C tines 
small  perturbation  parameter  e which  is  choser.  as  the  non- 
dimensional  thickness  of  the  plate  (set  equation  (3)  below). 
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^ 2.  Hooke's  Law  is  valid  between  stress  and  strain 

components . 

3.  The  coefficient  of  shear  viscosity,  at  the  median 
plane,  is  a constant  multiple  of  the  coefficient  of  normal 
viscosity  at  the  same  plane. 

The  first  assumption  is  similar  to  Goodier's  [8] 
assumption  in  treating  problems  of  beam  and  plate.  He 
assumes  that  the  components  of  the  stress  tensor  are  developable 
into  power  series  of  non-dimensional  thickness  parameter  e. 

In  treating  the  finite  deformation  problems,  this  method 
becomes  very  complicated  and  cumbersome.  It  seems  that  the 
same  purpose  can  be  accomplished  directly  by  dealing  with  the 
deformation  components.  We  thus  start  from  a point  at  which 
he  ends  up.  The  advantage  of  this  orocedure  is  that  it  avoids 
integration  of  some  systems  of  differential  equations  which 
may  be  quite  complicated  in  a finite  deformation  theory. 

Difficulty  arises  however  in  obtaining  desirable  orders  of 
magnitude  for  transverse  stress  components  for  the  isotropic 
material.  While  this  ooint  is  easily  explained  for  an 
orthotropic  material,  it  cannot  ue  explained  for  an  iso- 
trooic  material  without  considering  the  higher  order  terms  in 
the  expressions  of  the  deformation  ocmoonents . bucn  an  effort 
seams  to  be  hardly  worth  the  trouble. 

lhe  second  assumption  is  controversial  but  a commonly 
used  one. 

**  The  last  assumotion  need  not  be  included  into  the  theory. 

W ■ 

•is. 

The  whole  theory  can  be  carried  out  without  using  this  assump- 
tion. However,  unnecessaiy  coraoli cations  a^e  avoided  by  using  it. 
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The  present  theory  ends  after  obtaining  differential 
equations  of  the  lowest  order  pertinent  functions  in  the 
expansion  of  deformation  components. 

Five  partial  differential  equations  are  obtained  whose 
solution  is  believed  to  exolain  the  behavior  of  finite 
oscillations  of  visco-elastic  olates  to  a reasonable  degree 
of  approximation.  These  equations  are  non-linear,  and  contain 
all  previous  technical  plate  theories  as  special  cases. 

Boundary  and  initial  conditions  are  discussed.  It  is 
found  that  five  conditions  are  needed  for  every  part  of  the 
edge  surface  which  has  continuous  normals  possessing  first 
order  continuous  derivatives. 

Finally,  c.n  expression  for  the  strain  energy  is  given 
which  consists  of  two  parts,  namely:  Elastic  energy  and 

Energy  of  dissipation. 

2.  Components  of  Strain  Tensor 

•jo*' 

Let  x1,  x , x°  or  invariably  x,  y,  z be  the  rectangular 

cartesian  coordinates  of  a point  of  the  plate  before  de- 

1 2 

formation.  We  select  the  (x  , x ) -plane  as  the  median  plane 

3 

of  the  olate  and  x as  the  direction  perpendicular  to  this 

2 3 

"lane , so  that  x , x^,  xJ  makes  a right-handed  coordinate  system. 
3 

The  planes  x - ♦ h/2  will  be  taken  as  the  upper  and  lower 
faces  of  the  plate.  Components  of  the  deformation  vector  in 
rectangular  cartesian  coordinates  will  be  denoted  by  u1  or 
{ u , v,  w)  and  ccvariant  components  of  the  strain  tensor  by 
z • • or  (e  ...,e  }.  The  length  element  ds  in  the  deformed 

X J a A a y 

body  is  given  by  [9]: 
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, 2 

ds  * 

g^  dx^dx^  * ^ij  + 2eij^ 

(1) 

2ci  j 

- «»i  * ^4  (i , j - 1.2.3) 

<3xJ  dx1  dx1  dxJ 

(2) 

wherei 


1 when 
o when 


are  the  Kronecker  deltas , and  g^j  are  the  covariant  components 
of  the  metric  tensor  for  the  deformed  body. 

Throughout  the  paoer  the  usual  summation  convention  will 
be  used;  namely,  reoeated  indices  denote  summation  over  the 
range  (1,  2,  3).  We  note  that  there  is  no  need  to  distinguish 
covariant  mixed  and  contravar iant  corr.oonents  of  any  tensor  in 
cartesian  coordinates.  Hence  ■ eij,  u1  * u^. 

We  assume  that  u*tr^,x2,x3,t)  can  be  expanded  into  power 

A 

series  of  x°/h  * C: 


u/a-  e6[u0  ( ^ , T)  # t ) + t)  + e3C3u3  (£,,r|,  t)  ♦...] 

v/a*  e2[  vo(^,ri,t)  ♦ Cv-  (£,t],t)  + e3C3v3(|,n.t)  •‘...1 
w/a«  <wo(c,,r),t)  ♦ e2C2Wo  (c,.n,  t)  ♦ ( l .h  . t)  ♦ . . 


(3) 


where  rj,  C are  non-dimensional  coordinates,  a is  a typical 
length  and  e is  a small  oaraneter. 

x/a  - l,  , y/a  * rj  , z/a  = eC  , e«  h/a  , e «.  1 (4) 

The  explicit  forms  of  the  strain  coraioner.ts  are  obtained  by 
substituting  aquations  (3)  into  equations  (2).  Hence: 


i 
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e - e2[u  r + 4 w ,2  ♦ Cu-,  ♦ G(e4) 
xx  L o,£  2 o,\ 


cyy  ‘ ^'o.n  * k “o.r,2  ’ Cvl.^  * 0(e4) 


zz 


e2[^(Ui2  ♦ v^)  ♦ 2Cw2]  + 0(e4) 


e ■ 4e2[ u + v - + w -w  + C(u,  +v,  -)]  +0(e‘t) 
x.y  2 L o,r]  o,£  o,l  o,r)  l,n  1,£ 

cyz  ‘ 7 e(vl  * ''o.n1  * 0(e3) 

czx  ' 7 eiul  * vo.i)  *0(e3) 
where  indices  after  a comna  represent  dif ferentiation,i.e. 
uQ  ^ “ duQ/c)£  etc.  In  equation  (5)  the  second  order  terras 
represent  the  effect  of  large  riaf orraat ions  which  would  be  en- 
countered in  a second  order  theory.  The  second  order  terms 

occuring  in  the  expressions  of  e e , ®nd  e are  identical 

xx  yy  xy 

to  those  of  F&ppl-Karman-Timoshenko  theory  of  large  static 

deflection.  The  classical  olate  theory  assumes  u.  ■ -w  „ , 

1 o,t' 


> (5) 


v,  **  , hence  (c.  , z ) ■ 0(eJ).  Consequently  the 

i o,t)  yz'  zx 

transverse  shear  strains  e , , z . as  well  as  e , are  tal 

yz'  zx  zz 

as  aero.  The  present  theory  will  not  ignore  these  strain 


comoonents . 
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3 . Stress-Strain  Relations  for  an  Elastic  Medium  wi th  Internal 

Damping 


The  literature  contains  various  attempts  to  generalize 
the  classical  stress-strain  law  of  elasticity  so  as  to  include 
the  effect  of  internal  damping.  Early  discoveries  are  due  to 
0.  Meyer  [10].  W.  Voigt  [11]  and  K.  Sezawa  [12]  which  can,  in 
the  modern  notation,  be  expressed  as: 


ij 


- (\  + 0 


(^x  +2u'  ^)e ij 


0 - 


e 


l 

i 


(6) 


where  X.  and  p.  are  the  classical  Lame  constants  and  and  ' 

are  the  normal  and  shear  solid  viscosities.  Here  j are 

the  mixed  components  of  the  stress  tensor  which  are  referred 

to  the  curvilinear  coordinates  y*  of  the  deformed  state, 

while  strain  components  are  referred  to  the  cartesian  coordinates 

x*  of  the  undeformed  state. 

1 ? 

Coordinates  y and  y ' are  perpendicular  to  each  other  and 

3 

lie  in  the  tangent  plane  of  the  median  plane  and  y is  per- 
pendicular to  this  plane  and  makes  a right  hand  system  with 
y^"  and  y^. 

The  controversial  assumption  underlying  equations  (6) 
seems  to  be  as  good  as  any  one  in  use  relating  stress  components 
to  strain  components. 

Coefficients  are  related  to  Young's  modulus  E, 


Shear  modulus  G and  Poisson's  ratio  v by  the  following 


relations: 


\ « vE/(l-2v)(l*v)  . n * G (7) 

In  a similar  manner  E',  G'  and  v can  be  defined  by: 

- vE'  / (l-2v)  (1+v)  , p'  - G' 

Consequently,  if  one  defines : 

\ *K  ^ - H ♦ n'-JTf 

i - E * E'w  • g - ® * G w 

equations  (6)  take  the  same  form  as  those  of  the  classical 
elasticity  theory,  which  further  leads  to  a considerable 
simplification. 

We  note  that  the  strain  components  in  equations  (6) 
are  exoresned  in  cartesian  coordinates.  Hence  covariant, 
contravariant  and  mixed  comoonents  of  stress  tensor  need  not 
be  distinguished  for  the  present  approximation. 

When  a Voigt  solid  is  subject  to  oscillating  stresses,  the 
rate  of  dissipation  of  energy  is  proportional  to  the  square  of 
the  frequency  of  oscillation,  as  may  be  observed  from 
equation  (6)  for  the  one-dimensional  case.  Lord  Kelvin  [13] 
has  observed  that  the  rate  of  dissipation  of  energy  increases 
less  rapidly  than  the  square  of  the  frequency.  This  led  to 
the  generalization  of  Poynting  and  Thomson  [14]  and  more 
recently  to  that  of  Alfrey  [15,  16],  which  contains  all  other 
theories  given  before  as  a soecial  case. Alfrey's  stress- 

* L.  boltzmann  uses  the  superposition  principle  which 

seems  to  be  less  suitable  for  the  ouroose  of  generalization. 
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strain  relations  for  an  incompressible  material  are  expressed 
in  the  following  form: 


Poik  ‘ P ° Sik  * »eik  ' 0 ' T(oll  4 °22  4 °33> 

£ii  * C11  4 e22  4 e33  * 0 


m 

P - E 
i-o 


\ 


\ (10) 


J 


where  conatants  a^  and  are  to  be  determined  from  the 
experiment. 

We  now  generalize  this  equation  further  to  include  normal 
viscosities  in  the  case  of  compressible  material,  and  write: 


Poik  “ Re5ik  * 2Qeik 

l a1 

R "i-o  Ci  a?  ' 9 " Cii  " El1  +E22  +E33 


1 (ID 


j 


Here  P and  Q are  the  same  as  in  equation  (10) . We  now  solve 
°ik  ^rom  differential  equation  (11)  by  using  operational 
methods . 

Let  F(s)  be  the  Laplace  transform  of  F(t)  as  defined  by: 

F ( 8 ) - e"stF(t)dt  (12) 

Applying  the  Laplace  transform  to  the  expression  of  Po^  in 
equation  (11)  (i.e.  by  taking  Po^  for  F(t)  in  equation  (12)), 
we  obtain 

<Jik  * L(s)  6 (x1#  x2,  Xj,  s)  8^  +2M(s) 

^ij (*!*  x2*  x3 * s)  + N(x1#  x2,  x3,  s) 


(13) 
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where: 


L -i  m 4 n , in  • 

L(s)  - ( E c^s1)/(  E a, a1),  M(s)  - ( E b^s1)/^  a,s”) 
i*o  * i*o  i“o  imo 

N (s)  - (mE1D.si)/(  E a , s1)  , D,  - D. (x, ,x«,x«) 
i-o  1 i-o  1 1 11/3 

Taking  the  inverse  transform  of  equation  (13)  we  obtain 
aij  " 5i/o  Lo(t“T)*  6(xl'x2'x3'T)dr 


(14) 


+ Mo(t“T)  eij (x1.x2»x3'T)dT  * N0(xl'x2rX3't)  (15) 

If  the  system  has  no  initial  stress  then  o^.  must  vanish  with 
e^j#  thus  giving  NQ  « 0.  Furthermore,  writing  the  elastic 
parts  separately  in  equation  (15)  we  have: 

oi:j  - >.6bi.  + + Sij  £"  L(t-r)  6(x1,x2,x3,t)cIt 

+ 2 M(t-x)  (x^,x2,  x3,t  • dr  (16) 

where: 

L(t ) - L0(t)  -X  M(t)  - M0(t)  -n  (17) 

Thus,  we  arrived  at  V„  Volterra's  [17]  stress-strain 
relation  for  a medium  having  heredity.  Consequently,  the 
general  internal  damping  law  stated  by  equation  (11)  is  identical 
to  hereditary  damping  which  is  fr  iquently  treated  as  a different 
and  rather  more  general  phenomenon.  Equation  (16)  is,  of 
course,  identical  to  those  used  by  £.  Volterra  [16]  for  an 
elastic  medium  with  hereditary  damping. 
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A useful  form  of  the  stress-strain  relations  (16)  is 
obtained  if  we  use  the  convolution  multiplication  defined  by: 


F * i>  • Sq  F(t— t)$(t)  dr 


We  note  that  ? * 4>  * 4>  * F.  Equation  (16)  can  now  be 
written  si.aply: 


(IB  ) 


£e  5ij  * \ £ij 


(19) 


whexe  the  operators  \ and  q are  defined  by: 

* * 


X - X + L- 

ft- 


U + M* 


(20) 


Equation  (19)  has  exactly  the  same  form  of  elastic  stress- 
strain  relations  for  a homogeneous,  isotropic  media.  There- 
fore, substitution  of  X and  U for  X and  p.  of  classical 


elasticity  theory  in  all  equations  of  elasto-dynamics , gives 
the  corresponding  equations  for  a medium  with  internal  damping. 
Attention  must  be  given  to  the  fact,  however,  that  X and  p. 
have  time  derivatives. 

A more  compact  relation  would  be  obtained  if  we 
eliminated  6 in  equation  (19).  Writing  i ■ j we  obtain: 


0 “ (&  * T&)e 


(21) 


Subtracting  hydrostatic  pressure  a5^j  from  equation  (19)  , 
we  obtain: 


sij  “ * £ eij 


(22) 
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r 


O 


where  a and  e^  j are  the  deviatric  parts  -of  the  stress  and 
strain  tensors  and  are  defined  by: 


6 S 
T 5ij 


(23) 


sij  “ °ij  " o8ii'  *ij  “ eij 
It  is  important  to  note  that  the  functions  L and  M of  the  stress- 
strain  relations  (16)  may  be  discontinuous  functions.  Howevei , 
the  convolution  integrals  are  continuous  and  have  continuous 
derivatives  with  respect  to  time.  In  fact,  we  can  obtain  a 
general  expression  for  L(t)  and  M(t)  by  getting  the  inverse 
transforms  of  T and  IT.  In  the  most  general  case  we  may  have 
m ^ t , m<n,  which  corresponds  to  the  case  in  which  the  highest 
order  derivative  occuring  in  equation  (11)  is  that  of  o^.  For 
example,  the  Voight-Sezawa  type  of  damping  is  of  this  sort.  In 
this  case,  by  use  of  the  Dirac-Delta  function  and  its  derivatives, 

we  can  obtain  corresponding  X and  (i.  If  F(t)  a continuous 

* 41- 

function  with  continuous  derivative,  and  F(0)  » 0,  integrating 
by  parts  we  see  that  {*) : 

|| -Z1  8 (t-x)  ' F(x)dx  , 8 - (24) 


where  the  Dirac-Delta  function  5 (t)  is  defined  by: 
[oo  for  t ■ 0 


8 (t) 


■r 

lo 


and  X*  5 (t)dt  • 1 
for  t f 0 -cd 


(25) 


Hence,  equation  (6)  is  identical  to  equation  (19)  with: 


\ K \ ♦ X'8  (t)j 


u'§(t)^ 


(26) 


( 1 


If  the  difference  between  the  highest  order  derivatives  of  o 


ij 


Strictly  speaking,  this  equation  is  true  only  in  a limiting 


vt+e 


sense:  dF/dt  • lim  X 

e-*0  * 

See  also  L.  Schwartz  [19],  vol.  II.  p.  17. 


8 (t-x,e) f (x)dx  where  lim  8(t,e)  * S(t), 

e-eO 
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and  occuring  in  equation  (11)  is  k,  then  L and  M will 
contain  k**1  derivative  of  the  delta- function ^ ^ . We  -can,  in 
fact,  determine  and  discuss  all  oossible  forms  of  X.  and  a by 
discussing  the  roots  of  the  denominator  of  L and  M.  We  shall 
not,  however,  go  into  this  matter  here,  as  this  would  be  of  no 
interest  in  application  for  we  do  not  know  coefficients  a^,  b^ 
and  c^.  We  shall  be  satisfied,  then,  by  remarking  that  the 
commonly  used  expressions  of  L(t)  and  M(t)  are  of  the  following 

m t 

exoonential  type  [18]:  £ A.e  i . It  is  easily  seen  that 

i-o  1 

this  is  included  in  the  inversion  oA  L and  M.  In  fact,  if  the 
denominator  of  L and  M has  m real  distinct  root  and  (£,n)<  m we 
obtain  an  expression  of  this  sort.  Various  attempts  have  been 
made  to  determine  A^  and  in  this  expression  [20], 

It  may  be  of  further  interest  to  note  that  when  a^,  b^, 
and  c^  are  functions  of  space  coordinates  (x^,  Xg,  x^)  we 
obtain  a further  generalization  of  internal  damping,  leading 
to  L and  M of  an  inhomogeneous  media  which  are  also  functions 
of  coordinates  and  time. 

Equation  (19)  also  suggest  that  in  the  case  of  anisotropic 
alastic  media  we  can  use  extended  stress-strain  relations,  by 
simply  taking  the  asterisk,  »,  as  a subscript  in  all  elastic 
coefficients.  In  the  following  develooment  this  will  be 
the  procedure. 

^Mathematical  justification  of  the  use  of  the  8-function  and 
its  derivatives  may  be  found  in  3cnwartz's  book  [19]. 


4.  Plate  Stress-Strain  Relations 


We  consider  an  orthotropic  plate  having  a cylindrical 
symmetry  with  the  axis  of  cylinder  in  the  direction  of  the 
z-axis.  The  stress-strain  relations  can  be  written  in  the 
usual  technical  notation  as: 


c . xx  v x „ e . 1 

£xx  ~T~  r ayy  IT  °zz'  exy  TZ  °xy 

a v ^ , 

eyy  " “ E °xx  + “F  " S7  azz'  eyz  “ IB"  °yz 

Z z 

ezz  “ “ IT  °xx  “ F °yy  + ' exz  “ IB^  °xz 


*\ 

y (27) 


where  E,  G#  and  v are  the  Young's  modulus,  shear  modulus, 

and  Poisson's  ratic*  respectively,  in  any  direction  on  the 

x-y  plane,  E and  Gw  aie  those  for  the  z-direction.  v is 

the  Poisson's  ratio  corresponding  to  the  contraction  in  the 

x-direction  for  a unit  elongation  in  the  z-direction# 

Likewise  v is  the  Poisson's  ratio  reoresenting  the  contraction 
z 

in  z-direction  for  a unit  elongaticr.  in  the  x-direction. 

These  elastic  constants  are  related  to  each  other  by 
symmetry  relations,  since  the  strain  energy  must  be  a 
positive  definite  form.  Thus: 

v/E  - v/E,  G - E/2 (1  «■  v)  (28) 

X z z 

To  obtain  the  plate  stress-strain  relations  we  now  select  a 
particular  type  of  ortnotropy  which  is  consistent  with 


*1 

r 
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equations  (28)  : 

v » e2av  , E - e2aE  - e2E'  , G - e2yG  - «2G'  (29) 

X z z z 

Here,  a and  y are  positive  numbers,  and  v E and  G are 
independent  of  e. 

We  solve  components  of  stress  tensor  from  equation  (27) 
and  use  equations  (29) , hence: 


xx 


yy 


zz 


zx 


■ 7^2  (Exx  +veyy>  * 0(e4)-  °xy  * 2GExy  * 0(®4) 

* 7^  (eyy  ",cxx)  * 0(e4)-  °y2  “ 2e2(3'eyz  * 0(e5) 

’ rv’l'V***  * Eyy>  * <1‘V,EZ1  + 0(e6)' 

- 2e2  G #ezx  ♦ 0(e5) 


>(30) 


J 


We  thus  find  that  the  components  of  the  stress  tensor  are  of 
the  following  orders: 

(oxx-°y .y'V  ' 0(e2)-  ‘V0**’  ' °(e3)-  °m"  0le4)  (31) 

This  result  agrees  with  the  physical  situation  that  in  a 

shear  deformable  media  shear  modulus  G„  is  small,  and  in  a 

transverse  stress  deformable  media  E is  small.  The  olate 

z 

stress-strain  relations  (30)  can  be  so  adjusted  as  to  include 
a soecial  tyoe  of  internal  damping  by  replacing  E,  G,  and 

G'  by  E,  G and  Gf,  which  are  defined  by: 

* * * 

E = E ♦ e*  , G = G + g*  , G'  - G'  + g'*  (32) 

* * *- 

where  £,  G,  and  G'  are  the  usual  elastic  constants,  and  star(*) 
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multiplication  is  defined  by  equation  (18).  We  note  here, 
however,  that  in  this  particular  internal  damping,  \ is  a 


constant  multiole  of 


s- 


As  a matter  of  fact 


. _ 2v 

1 T=*r  $ 


, H 

*• 


2lf+v) 


G 

* 


(33) 


We  shall,  however,  keep  G'  as  an  independent  operator  through- 


out the  present  paper. 

The  simplification  introduced  by  equation  (33)  is  not 
necessary  for  the  analysis.  The  whole  analysis  can  be 
carried  out  without  this  assumption.  However,  for  the  sake  of 
simplicity  we  shall  use  the  assumption  contained  in  equations 

(33).  Thus,  we  have  two  functional  operators  E and  Q ' 

* * 


instead  of  the  three  representing  the  hereditary  character- 
istics of  the  medium.  We  note  that  many  important  practical 
classes  of  problems  are  contained  in  the  present  type  of 
internal  damping. 

In  the  present  study  we  shall  use  the  stress-strain 
relations  (19)  subject  to  (33),  which  are  valid  for  three- 
dimensional  eiastic  media  with  internal  damping.  We  shall, 
then,  indicate  the  terms  which  should  be  excluded  from  the 
remilting  expressions,  whenever  the  foregoing  type  of 
anisotropy  is  assumed,  in  order  to  obtain  the  formulas  which 
would  follow  if  the  plate  stress-strain  relations  were  used 
to  begin  with. 
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5.  Internal  Force  and  Moment  Corn-.  ->nents 


In  a plate  theory,  integrals  of  the  stress  components, 

as  well  as  their  moments  about  the  x and  y axes,  play  an 

important  role.  Thus,  we  define  normal  force  co  iponents 

(N  ,N  ,N  ) , plane  shear  force  component  N , trc  '.averse 
x’  y z xy' 

shear  components  (Qx,Qy) , bending  moment  components 

(Mx,My.Myx) , ana  body  force  and  moment  components  (Px, 

P ,P  ) and  (G  ,G  ,G  ) by  the  following  equations  (Fig.  1) : 
y z x y z 


(N  N N N ) ‘ 
1 x* y*  z*  xy' 

h/2 

" s 

-h/2 

h/2 

(Qx'  V “ *A\/2 

(o 

xz 

(Mx.  My.  Myx)  - 

rh/2 

-h/2 

^xv  " ^vx  ' 
* * 

Nxy 

xx*  yy'  zz'  xy 


xx'  yy'  yx 
•Nvx 


(a2e3)  IP  P P ) - J 'h/2  (F  . F . F ) dz 
y z -h/2  x y z 

(a3e4/12)  (G  G G ) - /h/2  (F  F F ) z dz 
x y z -h/2  x y z 


V (34) 


J 


where  F , F . F are  the  body  force  components  oer  unit 
x y z 

mass . 

Combining  equations  (5)  , (19)  and  (34)  v'9  obtain: 
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°xx  “ (Nx/h)  + (12/h2)  CMx  + 0(e4) 

Oyy  - (Ny/h)  ♦ (12/h2)  CMy  - 0 (e4) 

°zz  " ®2(q  + 2PC)  + 0(e4) 

°xy  * (Nxy/h)  + (12/h^  £Myx+  0(q4) 

°yz  " (Qy/h)  + °(e5) 

°xz  " (Qx/h)  + 0(e5) 

where: 

2P  “ [ozz(x,y,h/2.t)-  oz  (x,y,-h/2,t)]/e2-X.(u1  - 

* ' * 

* Tl,lj  * 2w2>  * * w2 


2q  “ Zozz(x,y,hl2tt)  ♦®w*(x,y,-h/2,t)]/e 


zz 


2X. 


Tu  r + v 
L o,l  o 


.r,  * $»o,i  * H.n2  * r,ui: 


+ V. 


2)] 


> 


J 


(35) 


(36) 


♦ (u-,2  + v,2  ) 


12  2 

If  we  eliminate  Wg  and  y(u^  + ) between  equations  (34) , 

after  using  equations  (19),  and  (36)  we  obtain: 

Vh  - TZ*  If  VS*  T “o.S2+v  (Vl*  K.^+T^T  1 

2 


N /h  “ — 9E[t  + iw  2+v  (uo  + iwrt  2)]- 

y l-v2*  °'^  2 °*T*  °#l  2 °»?> 


.-v 


Nz/h  - e q 


> (37) 


(u. 


+ V - + w^  - w ) 


Nxy/h  * TT&ni  '“o^  ’ To,S,  ‘ "o.S’o.t,' 

Qx/h  -«.V  (u1+  »3>l) 

Qy/h  -«eV  (v,  * w J 
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and: 


12Mx/b2  - 7^7  S <“1.5.  ♦ * -I#-  P 


2 2 
12M  /h2  - -2-,-  E U,  *»u,  ,)  * 2“S_  p 
y 2. -v 6 * ^ * h 1 » £ 2 -v 


> (38) 


12V/h' 


— - E (u,  + Y,  ?) 

2 (1+v ) * i'T5  1'^ 


We  note  that  the  terms  containing  p and  q in  equations  (37) 

and  (38)  appear  because  we  used  the  stress-strain  relations 

(19),  which  are  given  for  three-dimensional  isotropic  media. 

Instead,  if  we  use  the  plate  stress-strain  relations  (30),  we 

find  that  these  terms  disappear  from  the  expressions  of  N^, 

N , M . and  since  in  tthis  case  o = O(e^)  hence  p and  q 

y x y zz  • 

4 

become  of  order  e . See  also  [8].  In  this  case,  obviously 

Nz  ■ O(e^),  but  it  must  be  taken  into  account,  as  there  is  no 

lower  order  term  in  the  expression  of  Nz/h. 

The  introduction  of  a constant  k into  the  expressions  of 

Qx  and  Qy  is  due  to  the  fact  that  transverse  shear  stress 

distribution  across  the  thickness  is  not  constant.  It  is 

known, however,  that  e „ and  e„„  are  even  functions  of  z. 

yz  zx 

Hence  integrations  across  the  thickness  introduce  such  a 
constant#.  Such  a constant  was  first  introduced  by 
Mindlin  [4],  and  sucessfully  replaces  k ' of  Timoshenko's 
beam  theory,  [ 21]  and  5/6  of  the  Reissner's  plate  theory  [2], 
Mindlin  has  also  suggested  a formula  for  U in  his  paper. 
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It  may  be  of  interest  also  to  note  that  we  reached  the 

present  combination  of  E and  v in  the  first  terms  of  equation 

*• 

(37)  and  (38)  which  are  identical  to  those  of  the  classical 
theory,  not  by  using  two-dimensional  stress-strain  relations, 
but  by  using  the  three-dimensional  stress-strain  relations. 

6.  Equations  of  Motion 


The  equations  of  motion  of  an  elastic  body  referred  to 
the  coordinates  in  the  deformed  stage,  are  sufficient  to 
determine  the  exact  plate  equations  of  motion  in  a large 
deflection  theory.  Integration  across  the  thickness  in  the 
deformed  position  makes  it  possiole  to  introduce  membrane 
forces  and  the  plate  shear- force  and  bending  moment  components 
systematically.  In  the  present  theory,  the  plate  thickness 
is  assumed  to  remain  constant  throughout  the  motion. 

Equations  of  motion  in  curvilinear  coordinates  are  [9]: 


ilJ. 


pU1  - F1) 


(39) 


where  p is  the  mass  density  per  unit  volume,  F1  are  the 
components  of  body  forces  per  unit  mass  and  a^  are  the 
components  of  the  acceleration  vector.  Indices  after  the 
semicolon,  ; , represents  covariant  differentiation,  i.e: 


do 


ij 


U ^ 


rn 


i 


♦ n. 


,ia 


(40) 


Here  /T^  is  the  Christoff el  symbol  of  the  second  kind  and 
is  related  to  the  Christoffel  symbol  of  the  first  kind  ru.c 
and  the  components  of  the  fundamental  metric  tensor  g1^  and 


( 


t 
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r 


r 


^ij  *^e  c^e^orira<^  medium  by: 


rk.rk.gtV 

ji  ij  ij ,a  dxJ  dx1  dxa 


(41) 


dxJ  dx 

TJslng  equations  (1),  (5)  and  (41)  we  find  that  the  components 


°f  qj  k are: 


n.i  " (e2/2a)  tu°.t  * T"  0.52*  Cui.i).r 


r " (oZ/2a)  (v  + w w ♦ Cv-,  rx  ) , 

^■2  °*U  °.n  i#U 


r - (e/2a) w ..  , r " r 

11.3  °'“  12.1  21,1 


(e^/2a) (u 


o.i 


* iwo.i2  + 


r - r - U2/ 2a)  (T  * i«o  a2  " Cvl  n*  f 
12,2  21,2  °'^  ^ 1,11 


r - r • (e/ 2a) w . , r - r 


12,3  21.3 


O.lT,' 


13,1  31 


- (e/2a)u, 

.1 


2 1 «_ \ / 1..  2 


r -r  “ (e/2a)v,  *,  r ~r  ■ (e^aM-iu^ 

13,2  31,2  A#t  13,3  31,3  * 1 


- - 2;w2) 


i(42) 


r 

22,1 

- (•*/*•) (»0<n. 

+ w rw  + 

0,1  o,T|T) 

Cu,  ). 

1 ,nh 

r 

22,2 

- <•*/*»> <*0., 

+ Jlw  ^ + Ty 

* O.T)  Vl, 

n),-n  ' 

r 

■ (e/2a) v , 

r - r - 

(e  / 2a)  u. 

9 

22,3 

23,1  32.1 

1 ,n 

r 

23,2 

- l~  -(e/2a)v1 
32,2  1 

„ , /”  -r  - 

'*>  23,3  32,3 

(e2/2a)  l^2 

4vi 

♦2Cw„)  ...  /'  — 

r , r -r 

, r ■ (e/a) w9 

* 33.1 

13,3  33 , 2 23,3 

33,3 

“l 


V: 


“I 
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The  terms  of  the  lowest  order  in  o1^  . are: 

9 J 


lj 

_ da^ 

tj2j  - 

do2-* 

; j 

dx-' 

' ;j 

dxj  ' 

3 j 

. *°3i 
> 

♦ r3  a11  ♦ 

2T3  o12  + 

r3  a22 

» j 

dx*1 

11 

12 

22 

3 3 

Thus,  in  the  present  theory,  we  only  need  to  retain  f~  , f~ 

11  12 

and  f~,  A more  general  theory  will  undoubtedly  contain  more 
22 

of  these  coefficients.  Therefore,  equations  (42)  may  be  useful 
for  that  purpose. 

Using  equations  (1),  (41)  and  (42)  we  find  that: 


♦ 0(eJ) 


l~^m  r + 0(e3)  - (e/2a)  w __  + 0(e3) 

11  11,3 

r3  - r + G(e3)  - (e/2a) w + 0(e3)  (44 

12  12.3 

r3  “ r + 0 (e5 ) = (e/2*)w„  + 0(e3)  j 

22  22,3  °'T1T1 

Components  a1  of  the  acceleration  vector  in  the  cartesian 
coordinates  x1  are  given  by: 

? i 

_i  _ on 


Hence  by  equations  (3)and  (45): 

a1  - ax  * a e2(uQ  + Cu^)  + C(e4) 

a^  * ay  * a e^(vQ  + Cv^)  + C(e  j 

3 **  3 

a - a m a e w 0 (e  ) 
z o 


Microscopic  equations  -f  met  on  (39)  t:  us  become: 
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do 


xx 


da 


3x 


Ty*- 


d a 

xz 

<3z 


pux  - 


Fx> 


'V 


ox 


da  da 
TL  + _ja 
dy  dz 


>(47) 


da  da  da 
zx  zy  . zz 

dx  dy  dz 


r3o  ^ 
11  xx 


2r3o  ♦ r3o 

12  xy  22  yy 


P(az‘Fz} 


Mass  density  p is  given  bv 

p - pc(l  - 20  + 462  - 8 93)1/2  - pa[l  - 0(e2)  ] (48) 

where  0 , 02,  ®3  are  *he  strain  invariants  of  which  0 *e^ 

2 

•0(e)  is  the  lowest  order.  Hence  in  the  present  theory 
p * p , that  is  the  change  in  mass  density  during  the  de- 
formation is  negligible.  The  force  and  the  moment  equilibrium 
of  a plate  element  with  thickness  h,  cut  off  from  the  plate 
by  drawing  normals  to  the  deformed  median  plane,  is  commonly 
used  [5]  to  reduce  the  microscopic  equilibrium  conditions  to 
macroscopic  equilibrium  equations.  This,  of  course,  represents 
an  approximation  which  replaces  the  correct  equilbrium  condition, 
namely,  that  moments  of  all  order  across  the  thickness 
should  be  zero.  The  foregoing  practice,  however,  only 
insures  zero  moments  of  the  first  and  zero  orders.  This 
common  practice  is,  however,  in  accord  with  the  present 
tneory.  Hence  integrating  equations  (47)  across  the  thickness, 
we  obtain: 


dN  dN  h/2 

W~  * [oxz-h/2 


h/2 

4/2  P(ax 


Fx)dz 


►(49) 
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dH 

T? 


aN 

♦ — L 

~>y 


[ay2] 


h/2 

-h/2 


h/2 

/-K/2PUV  ’ V* 


I 

I 


aQ„  ag  h/2 

“ST  + “3y  + ^zz1  + 

-h/2 


11 


r3  nx  ♦ 2 r n ♦ r\ 

*.  -,n  xy  22  y 


12 


h/2 

/ P(a  - F )dz 
-h/2  2 2 


y (49) 

J 


Next  we  multiply  equations  (47)  by  z and  integrate  across  the 
thickness.  Hence: 


air.  aNi  . 

-Ji!  * ?°xz(*.ir'h/2.t)  + 0xz(x.y.-h/2,tJJ-Qx 

h/2 

* / p (^  - Fx)  z dz 
-h/2 


aM„_  aM 


dr  + •3yZ+  7C°yz(x,y,h/2,t)  ♦ oyz(x,y.-h/2,t)]-Qy 
h/2 

" ^ tfav  " FJ  z dz 

-h/2  y y 


* 


h/: 


h/2 


(50) 


-h/2 


°xz  2 


dz)  + 17  ^Llfl0yz  2 d2)  * ?Ozz(*'y»h/2,t) 


-h/2 
h/2 


♦ °zz{x'  y>h/2,t)]-Nz  “/  P<az  * Fz)z  dz 

-h/2 


We  assume  that  the  surface  shear  is  zero: 

°xz'X,y'"h/  2 * 1 } " ay2(x.y^h/2.t)  = 0 (51} 

Using  equations  (37),  (38),  (44),  and  (46),  in  equations 
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r 


(49)  and  (50)  we  obtain: 
.3 


7-E*Cuo,u’vv°.’ll  + 7 (wo.i2<',uo.T,2).ll 

3 

♦ T“ — «3q  r ♦hh;  > £ [u  ♦ V r + (w  _W  ) 1 

1-V  ,1  2 U+V  ) *>•  O.TJT]  0,&T)  0,1  0,T),T)J 


A 


♦ 0 (e5)  - pa2c3(u0  - Px)  ♦ 0(e5) 


3 

+ ',uo,5n  + 7(wo.n  1-n^  * Tr\T  ® J 

* 7n*v)EJv^  * ’o.u  * (wo.:wo^U]  * 0(e5) 


Da2e3(‘v0  - P ) ♦ 0(e5) 


* ’I.n  + wo.U  * %,o.nn)  * ^ 

+{r^  + ^o.t2  *v(vo.r, + K,n2)j 

4 1 

♦ rv  .-e^q  [ w __  + r-rTT — r E (u  ♦ v - ♦ w -w  ) ] 
1-v  J o,ll  L U+V  ) *V  OrTl  o,^  0,1  o,T]  J 

w _ *( EJv  ♦ i w 2 * V (u  - ♦ i w v2)] 

°»ST)  1,1 -v  2 **"  2 0,T]  o,\  Z o,\ 

lq  "I  w ♦ 0(e6)  * oa2e2w  + 0(e4) 

J O.IT]  O 


>(52) 


, v 4 
+ T=T  e 


12 ( 1— v ) 
(u 


3 3 3 

~77  E*(ul,U  + vvl,lV  + ffTT-vT  p,l  + 24(l*v)  ' E« 


\ 


v,  _ ) "#i3G  (u,  •*•  w _)  ♦ 0(e5) 

l,fcTj  * 1 o»l 


pa2e3  (*u^  - Gx)  ♦ 0(e6)  , 
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7TE*  lvl.nn  * Vul,ltl)  4 TO77)  p, 


12(l-v2)  * 

3 


* 2j(*+v)  E.  <ul,n^  4 Tl.u’  4 wo,») 

♦ O(e^)  “ P ~ Gy.)  ♦ O(e^)  , 

2T(T+v T E J ( 2 w2wo , 1 ^ + (2w2wo,t)^  * 0 ' 


V (53) 


Equations  r52)  and  (53)  are  the  equations  of  motion. 

We  note  that  to  obtain  these  equations  we  used  three- 

dimensional  stress-strain  relations.  These  equations  can  be 

2 2 

reduced  further  if  we  take  e q and  e p for  q and  p respectively, 

thus  leading  to  a theory  in  which  two-dimensional  plate  stress- 

strain  relations  are  employed.  In  this  case,,  all  terras  in- 

2 

volving  q and  p drop  out  except  the  2e  p term  of  the  last  of 
equation  (53) . We  also  note  that  the  inertia  term  in  the 
latter  equation  is  of  0(e  ).  This  point  will  be  cleared 
when  we  used  perturbation  in  time  t later.  As  the  last  of 
equation  (51)  is  of  O(e^),  we  ignore  this  equation  in  the 
present  theory. 

7.  Reduction  of  the  Equations  of  Mot i on 

The  third  equation  of  (52)  may  be  written  in  the 
following  form: 


1 
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*e2G'  (u,  - + ▼,  J « -*e2G' Aw  + pa2w*  - 2p  - e2N, 

* 1*T>  * ° ° 

2 2 

•2hN  - N w + 2N  w _ + H w , A “ — — ir  ♦ -i_ 

x xy  y o*7!7!  ^^2  dt] 

From  the  first  and  the  second  of  equations  (53)  v^  and  u^ 
are  eliminated  with  the  use  of  equation  (54) . Hence: 

e aUi  - kg-  ux  - ^ Pa2(;\  - v ii| 

* 24*1— r 5 <P  ‘*®'Awo,i  "V 


(54) 


♦ (2ve*#G#-  E)  ep  ,/12(l-v)  - 0 

* * 


>(55) 


Avl^G#Yl‘  IlPa2(vl  ■Gy)'^£#Wo,^ 


* 24*l-v J l C(P  a2/*2)wo,T1  ~ “V 

♦ [ 2ve2  AC  G'-E  ] ep  /12 (1-v ) - 0 

» * ,v 

where  we  used  the  following  properties  of  the  operators 
E and  G ’ : 


J 


E G1  ■ 8'E  , A [E  #(*,!,, t>]  - E |A  . 

******  * ^ 

^ [E  *(l,»l,t)]  « E 


> (56) 


which  can  be  proven  easily  from  the  defining  equations  (32). 
We  can  also,  by  part  integration,  prove  that: 
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(E  i>)  , - E (* . ) - e (t)^(l.r).O) 

*’  ft  9 ^ 

(E  *)  - E(^  ++)  - e(t)*(£,r],0)  ♦ e(t)i(^,r)#o) 

* ' * ' 


>(57) 


These  identities  will  be  used  in  the  following  analysis. 

Differentiating  the  first  of  equations  (55)  with 
respect  to  £ and  the  second  with  respect  to  rj,  adding  the 
results  to  each  other  and  using  equations  (54),  (56)  and  (57) 
we  obtain: 


e2G'D 

* * 


A2w_ 


pa2  (D 

* 


^tfe2G')Aw^  - pa4h 


tfe2G' 


w 


- P13  • w,-  ♦ 2[ -D  + K.2V]Ap  42h3KG'p 

4 P-h-Ea2  P-  «*D  AN  » h3  *.V  N ♦ V 

D * ft  U 

4 eT!^  « *2G- (Gx.*  4 4*#tY. 


>(58) 


♦ u 


i.l  * ’l.n1  "° 


J 


where : 


® “ 12(l-vz)  * 

[Y,  - Y (t) ^(^,r),o)  ♦ Y(t)i(^,T),o) 


(59) 


Here  the  operator  D takes  the  place  of  the  bending  rigidity 

* 

of  the  plate  theory  without  internal  damping.  Cy.»«0  is 
known  when  the  initial  conditions  are  given.  Equations  (54), 
(55)  and  (58)  are  the  equations  of  motion  for  rotatory  and 
transverse  motion.  In  case  of  zero  internal  damping,  [y,...}20 
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and  star,  indexed  quantities  G',  E and  D must  be  replaced 


by  G',  E and  D.  Hence,  we  can  divide  these  equations  through- 
2 

out  by  Af  e G*  leading  to  plate  equations  with  no  internal 
damping.  This  is  done  only  for  equation  (58)  leading  to: 


-DA2w  ♦ , Dga!  j A"  . 4h  - 

° 12  tetQ’  ° ° 


-(pW/12^«2G')‘v  ♦ 2[ 2 — + — ]Ap 

° K»ZG'  1 2 ( 1-v  ) 

+ (2h3/e2)p  + (ph3a2/12*e2G')  2p  - (D/*G')AN 

+ h3N  + (ph5/ 12  K e2G')N  + (pa2h3/12) ^ - ♦ Gv 

x 0 s y 9 m ' ^ 

D - Eh3/ 12 (1-v 2) 


>(60) 


Equation  (60)  takes  the  place  of  the  claesical  plate  equation: 


-D  A2w_  - oa4h  w*  + 2a2h  p ♦ h3N  = 0 (61) 

which  are  four  terms  of  equation  (60) . The  terms  containing 
G*  represent  the  effect  of  shear  deformation  and  the  terms 
which  are  differentiated  with  respect  to  time  and  which  con- 
tain  (ph  / 1 2 ) as  a coefficient  represent  the  effect  of  rotatory 
inertia.  The  effect  of  the  large  deformation  is  contained  in 
the  expressions  of  N,  since  N , N , and  N , given  by  equation 

a y j wy 

(37),  are  multiplied  by  various  derivatives  of  w in  the  expression 

(54)  of  N leading  to  a non-linear  equation.  Body  moment 

components  are  taken  into  account  with  terms  G„  , and  G.  - 

* t \ y.ri 

The  term  [ (vh  /12(l-v)]Ap  is  rot  encountered  in  the  classical 
theory.  This  is  due  to  the  transverse  stress  component  which 
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r* 


is  ignored  in  the  cla:  sical  plate  theory.  The  foregoing 
term  is  the  moment  of  the  plane  forces  due  to  contraction  in 
the  x and  y direction  in  the  presence  of  surface  pressure 
p.  Hence:  In  order  .o  obtain  classical  plate  equation . we 


must  let  G*  — »oo.  ph  / 1 2 
and  p • 


_o , while  pa^h  “ finite  valu< 


JLJL' 


Axial  disolacements  u and  v satisfy  the  first  and  second 
• o o J 

of  equations  (52)  which  may  be  transformed  into  simpler  forms: 


E tsx>  - 2(l+v ) oa2to  * E (w  Aw  . - w _ Aw  ) 

*•  ° ° o.h  o,l  o.l  0,T) 


♦ [-  P * P ■} 


(62) 


E A©  - pa2  ( 1-v 2)  6 * -E  [ w ,Aw  _ ♦ w Aw  ♦ (Aw  ) 2 

^ ° ° ^ °.T)  o#1!  o 

♦ (l+v)(w  - 2 - w T?w  )]  - v ( 1+v ) Aq  ♦ pa2(l-v2) 


- Py^ 

where  go  and  ©^  are  the  rotation  and  the  dilatation  defined 
o o 

by: 


(63) 


6 * u « ♦ v 

O 0,1  0,7)  • 


(0  “ u - V - 

3 °,*1  0,1 


(64) 


Between  u , v , and  6 , to  dual  relations  exist  which  are 
o o o'  o 

expressed  by  equations  (64)  and  (65); 

AU.  “ 0 r + CO  ^ AV  “ 0 “ <0  t (65) 

o 0,1  0,T)  o 0,T)  0,1, 

We  note  that  when  the  right  sides  of  equations  (62)  and  (63) 

are  taken  as  zero  and  E is  replaced  by  E,  we  obtain  the 

plane  rotational  and  dilational  wave  equations  of  the  small 

deflection  theory.  The  right  hand  sides  of  these  terms 

represent  the  effects  of  finite  deformation,  body  forces  and 


rv 


-31- 


n 


n 


the  effect  of  plane  contraction  due  to  a load  q. 

In  application,  boundary  and  initial  conditions  are  usually 

given  i cerms  of  planar  displacements  uQ  and  v , hence  it  is 

useful  to  express  differential  equations  in  terms  of  u and  v , 

o o 

This  is  done  by  first  expressing  the  first  two  of  equations  (52) 

in  terms  of  u , and  v , <a  and  then  eliminating  o>  with 

o'  o , tj  o'-  0,5,'  ^ o 

the  use  of  equations  (62),  Hence; 


2 2, 


A 'A  uQ  + v(l+v)A'  q ^ ♦ (1+v ) pa  [e,o>Q  ] * A 
Q + v (1+v ) A'  q - (14v)2pa2[i#«o  - B* 


A 'A*  v 


) 


(66) 


where: 


A*  f 


E^Af  - pa  2 ( 1-v  2)  f ' 


A'  f « E^f  - 2paz ( 1+v ) f 
A*  - -Ia' 


^'E*C(“o.l2  * v*0.Al  * (1'v)(wo..lwo.n).^ 


> (67) 


W-o.n^o.t.  ' wo,iAwo,ii''.i"f’a2*1”',2*A,p3 

* * py.^ 

B*  is  obtained  from  A*  by  interchanging  x,£  and  y,T],  respectively. 

Equations  (66)  are  the  equations  of  extensional  vibrations 
of  visco-elastic  plates  including  the  effect  of  finite  trans- 
verse oscillations.  It  is  interesting  to  note  that  the  terms 
containing  transverse  deflection  wQ  are  completely  separated 
from  the  axial  displacements  uQ  and  vQ.  Thus  no  coupling 
effect  is  intvoduced.  The  coupling  effect  between  uQ  and 
v is  introduced  only  with  the  brackets,  Te,  'o  1 and 
re,  o)  rl  which  are  the  results  of  initial  conditions. 

Moreover,  this  coupling  effect  is  introduced  with 
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the  rotation  u>0.  Those  terms  disappear  for  purely  elastic 
media  where  memory  function  e 2 o. 

Equations  (66)  are  linear  ir.  uq  and  v „ Consequently, 
once  we  obtain  wQ,  the  right  hand  members  A*  and  B*  will 
be  known.  Therefore,  the  extensionai  Vibration  problem  will 
be  reduced  to  solving  the  inhomogeneous  linear  partial 
differential  equations  (66)  . 

For  elastic  plates  with  zero  internal  damping  equations 
(66)  can  be  written  asr 


A2u  - (1+v)  (3-v)E‘1pa2Au’  + 2(l+v)  (l-v2)E“2p2a4’u’ 

o o r o 


+ v (1+v  ) E_1Aq  « - 2v  (l+v)2E_1pa2q  , 

9 b 9 b 


A E 


-2 


> (68) 


A2vq  - (1+v)  (3-v)E‘1pa2Av’o  + 2 (1+v  ) (1-v  2)  E'2p2a4’*v” 

+ v(l+v)E'1Aq  - 2v  (l+v)2E"1pa2q  - B E~2 

9 M #*1 

where  A and  B are  obtained  from  A*  and  B*  by  taking  E for  E^. 
Equations  (68)  are  the  generalizations  of  the  equations  of 
extensionai  vibrations  of  the  classical  theory.  In  the 
classical  theory  terms  containing  A,  B,  and  q are  ignored. 


8.  Perturbat ion  in  Time 

The  above  differential  equations  are  extremely  complicated 
and  do  not  lend  themselves  easily  to  integration.  Simplification 
can  be  made,  however,  if  some  additional  assumptions  are  made. 

For  examole,  in  the  application,  plates  have  finite  boundaries. 
The  stress  or  deformation  wave  in  the  plane  of  the  plate  has 
a very  large  wave  velocity,  Consequently,  a few  milliseconds 
after  the  motion  begins,  the  stress  wave  will  be  felt  every- 
where on  'che  plate.  Thus,  in  problems  where  the  transient 
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state  in  the  planar  wave  propagation  is  unimportant,  a great 
deal  of  simplification  can  be  effected  by  making  perturbation 
in  time.  Also  assume  that  the  body  forces  and  moments  are 

4 

of  0(e).  Hence  let: 

t = et  (69) 

The  first  two  of  equation  (49)  can  then  be  reduced  to: 


N*.l  ’ N*y.i 


* Ny.n  ' 0 


(70) 


since  all  terms  in  equation  (49)  are  of  order  e , while 
inertia  terms  are  of  order  e^  , (See  equations  (52)). 
Equations  (70)  are  satisfied  by  the  Airy  stress  function 
P0U.T)ft): 

N 


e2F 


o.hh 


N*y  * ’Ad,  ' 


"y  ■ e2po.U 


uQ  and  vQ  contained  in  the  expressions  of  Nx,  Ny  and  Nxy 
given  by  equations (37) , can  be  eliminated,  thus  leading  to 
compatibility  equations: 

(Ny  - vN)  ^ ♦ (N  - vN  ) „ - 2 ( 1+v ) N - - 0 

Substitution  or  equation  (71)  into  (72)  gives: 


AZF  * vhA  q ♦ hE_(w 


- w r T w ) 

O.U 


(71) 


(72) 


(73) 


-*-•0,1^ 

The  equation  for  the  transverse  vibration  is  obtained  by 
combining  equations  (58)  and  (71).  The  only  change  in  equation 
(58)  will  be  in  the  expression  of  N,  which  takes  the  form: 


hN  m F w * r - 2F  - w » + 

o,r\r)  o,H  o.It] 

Of  course,  we  also  set  G =.  G - 0, 

^ y 

containing  derivatives  with  respect 
respect  to  t,  by  using  equation  (69) 


F ..w 
o.ll  0,T1T) 


(74) 


and  transform  the  terms 
to  time  to  derivatives  with 
. Equation  (58)  with  N 


given  by  (74)  and  equation  (73),  are  solved  simultaneously  to 

determine  F and  w » Once  this  is  done  we  can  then  use  the 
o o 

first  two  of  equations  ( 37 )or  (66)  to  determine  uQ  and  vq  since  N^, 

N and  N will  be  known  by  equation  (71) „ Deformation 

y xy 

components  u^  and  v^  can  now  be  solved  from  equation  (55) 

We  note  that  equation  (73),  without  the  term  vhAq, 
and  equation  (58  with  all  terms  ana  the  internal  damping  taken 
as  zero,  except  the  first  and  the  ninth  terms,  gives  the 
large  static  deflection  equation  given  by  the  Foppl-Karman- 
Timoshenko  theory  (») , Thus  the  present  theory  generalizes 
this  theory  to  dynamic  problems.  Moreover,  it  extends  the 
theory  to  take  into  account  the  shear  deformation  rotatory 
inertia,  inertial  external  loading  on  the  surfaces  and  the 
boundaries  of  the  plates  for  visco-elastic  plates.  Various 
special  cases  given  below  are  suggestive  in  simplifying  these 
equations  further  for  practical  problems, 

a.  Zero  internal  damping 

In  this  case  equation  (60  becomes: 

-DA2wq  ♦ (pj^  * 5§^wo,tt  " P*2h3vfo,TT 
-(pVV/12«G')w0  TTTT  * 2(-^__  ♦ yv^^Ap  *2ha2p 

+(ph3a2/12*G')2p  _ - (D/#G') AN  ♦ h3N  ♦ (ph5/12K3')N  -C 

9 **  9 

where  n is  given  by  equation  (74). 

bo  Gmail  deflection  and  ( a) 

If  we  use  equation  (75),  and  in  all  the  equations  which 
follow  take  N as  giver,  by  the  second  of  equations  ( 54)  , where 

(*)  The  theory  was  given  first  by  Ftfppl  [1], 
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N , and  N are  constants,  we  obtain  small  deflection  theory, 
xy  y 

c.  Rotatory  inertia  omitted  and  (a) 

3 

This  means  dropping  all  terms  containing  ph  / 1 2 in 
equation  (75) , hence: 


-DA2wo  ♦ (Dpa2/*G')Avo>TT  - pa2h3-Q  TT 
2 3 

+ 2(— + 12HT-V AP  + 2hfl2P  ‘ (D/^GOAN  + h3N  - 0 (76) 

Rh  G ' 


If  we  further  put  ^ ■ 0 to  obtain  static  plate  deflection,  we 
obtain  Reissner's  equation  (70)  given  in  [22],  except  for 
the  term  vh  /12(l-v).  As  explained  before,  this  terra  does  not 
appear  if  one  starts  with  the  plate  stress-strain  relations. 


d.  Shear  deformation  omitted  and  (a) 

By  lettinaAfG'  — >ao  in  equation  (75),  we  obtain  the  plate 
equation  in  which  shear  deformation  is  negligible: 

•“2-o  ♦ (ph5/12)AwO  TT  - pa2h3wO  TT  ♦ 5$^  AP 

+ 2ha2p  h3N  - 0 (77) 

Further,  if  we  set  ph"/I2  * 0 we  obtain  the  plate  equations  cf 
Foppl-Karmcfn-Timoshenko  theory  except  for  the  term  containing 
vh3/12(l-v)  whose  presence  is  explained  above. 

e,  Small  deflection  with  axial  stresses  omitted 

By  setting  N*0,  equation  (76)  reduces  to  the  equation  given 
by  Uf lyand-Mindlin  theory  [3.  4]  except  for  the  term  containing 
vh3/12(l-v)  . 


\ 
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9.  Boundary  and  Initial  Conditions 


Let  nj  be  the  exterior  unit  normal  to  the  edge  surface  E 
and  be  the  stress  vectors  acting  on  E.  Then  the  boundary 
conditions  on  the  surface  tractions  consist  of  giving  on  E 
at  all  times: 

T*  ■ prescibed  on  E (78) 

For  the  present  theory,  there  is  no  distinction  between  contra- 
variant  stress  vectors  T*  and  covariant  stress  vectors  T^o^n-^. 
Let  s be  the  tangent  to  the  boundary  curve,  C,  which  is  the 
intersection  of  the  deflected  medium  plane  with  E.  Then  n, 
s,  z make  a right  hand  system c Let  0 be  the  angle  between  n 
and  x ard  ^ - 0 the  angle  between  n and  y,  Fig.  2 . We  first 
transform  0*-*  components  to  the  n,  sf  z system.  Afterwards, 
we  obtain  the  plate  boundary  conditions  on  the  membrane  forces 
by  integrating  T1  across  the  thickness  from  -h/2  to  +h/2. 
Boundary  conditions  on  the  bending  moments  are  obtained  by 
multiplying  T^  with  z and  integrating  with  respect  to  z 
across  the  thickness.  Hence: 

Nn  * Y(Nx  + Ny}  + 7(Nx  ' Ny)  cos  2 6 + N sin  26 


N 


ns 


i(N  - N )sin  20  + N cos  20 

4 y a aj 


(79) 


Q__  - Q,  cos  0 ♦ Q sin  0 

ii  £*  a y 

IL  • i(M„  ♦ M ) + i(M,  - M ) cos  20  - M sin  2© 
n z x y i x y xy 


M 


sn 


X - Vin  20  + yxy  cos 


20 


(80) 


We  note  that  like  in  Reissr.er's  theory,  [2]  we  have  three 
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boundary  conditions  (80)  in  case  of  zero  membrane  forces.  In 

the  classical  plate  theory, of  course,  the  condition  on  Q is 

n z 

absent.  In  the  presence  of  membrane  forces,  altogether  five 

boundary  conditions  must  be  satisfied. 

In  case  the  surface  tractions  are  unknown  Lut  some  support 

conditions  are  given,  depending  on  the  type  of  support,  we 

need  to  prescribe  deformations  uQ,  vq,  wq,  u^,  and  their 

derivatives  in  various  directions.  These  conditions  are  easy 

to  express,  from  the  meaning  of  the  deformation  components.  Of 

course,  zero  derivatives  of  these  functions  in  any  direction 

along  the  boundary  express  the  clamp  conditions. 

In  general,  we  might  have  a mixed  condition  at  the  edge 

involving  some  conditions  on  traction  and  some  conditions  on 

deformations  and  their  directional  derivatives.  Altogether, 

five  independent  mixed  conditions  are  needed  at  a oart  of  boundary 

on  which  the  normal  n is  continuous  and  has  first  order  derivatives. 

Initial  conditions  are  obtained  by  orescribing  the  dependent 

variables  and  their  time  derivatives  of  first  order. 

An  exhaustive  study  of  the  boundary  and  initial  conditions 

can  oe  obtained  by  using  the  variational  orinciple;  which, 

however,  seems  not  worth  the  trouble. 

Finally,  we  give  the  expressions  of  M , M ,srd  M . 

x y xy 

Hence  M and  M can  be  obtained  in  terms  of  w^  by  means  of 
n sn  o 

equations  (80),  sc  that  the  differential  equation  for  w_.  can 
be  treated  without  reference  to  the  other  deformation  components. 
This  can  be  done  by  combining  the  last  two  of  equations  (37)  with 
equations  (38).  Hence: 
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If  the  plate  stress-strain  relations  were  used,  the  terms  con- 
taining p would  drop  out.  In  case  of  zero  internal  damping  one 
can  also  cancel  the  coefficients  /{hG_ ' . 


Ajh  - -D(w  „ «■  v„  i + (a/h) 3 (D/AfhG' ) (Q  ,*vQ  ) 
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In  these  expressions  the  terms  containing  wQ  are  identical  with 
those  of  the  classical  plate  theory,  while  the  terms  containing 
Qx  and  are  due  to  shear  deformation.  Equations  (82)  are 
the  analogues  of  Reissner's  equation  (10)  of  [2], 
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10 . Strain  Energy 


The  present  form  of  the  stress-strain  relations  (equations 
19)  implies  that  the  expression  of  strain  energy  V is  of  the 
following  form: 

V « / Udv  (83) 

(v) 

where  U is  the  strain  energy  density  per  unit  volume. 

(84) 


U 


1 _ij_ 
7 0 e-’ 


ij 


This  form  of  the  strain  energy,  with  stress-strain  relations 
given  by  equations  (]9  ) , is  identical  to  that  of  the  infinitesimal 
theory  in  the  case  of  zero  internal  damping.  This  form  is  of 
course  valid  for  infinitesimal  strain  but  large  deformation  and 
rotations  [23],  which  is  the  basis  of  the  present  theory. 

Combining  equations  (19),  (20),  and  (84),  we  obtain: 
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U ♦ U . 
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Ue  " + 2n>0i 
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(85) 


Here  U0  is  the  elastic  energy  density  and  U^,  the  dissipative 
energy  density.  Quantit  es  9 and  0£  are  the  first  and  second 
strain  invariants,  of  which  the  first  is  given  by  the  second  of 
equations  (6)  and  the  second  is  given  as: 

► 1 O 

(86) 
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For  the  Voight-Sezawa  type  or  internal  damping  they  are 
given  by  the  second  members  of  equations  (26)  namely: 


L * \'£  . M - n'B  . (87) 

9 ^ 9 ^ 

We  now  use  the  expression  (5)  of  e . ; in  the  expression  (85) 

* J 

of  U and  integrate  with  respect  to  C across  the  plate  thickness, 
then  use  equations  (36)  to  (38).  The  result  is: 
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Here  and  are  the  strain  energy  densities  per  unit  area  due 

to  the  membrane  force  and  bending  moment  components,  respectively. 

Subscript  e at  the  end  of  a bracket  indicates  the  fact  that  only 

the  elastic  ->art  of  the  expression  inside  of  the  bracket  must  be 

used.  We  note  that  0^  and  Ujt,  consist  of  two  parts,  namely,  the 

elastic  energy  and  dissipative  energy  due  to  damping  of  which  the 

first  is  obtained  by  considering  the  elastic  parts  of  the  membrane 

force  and  bending  moment  components  wh’ch  are  multiplied  by  the 

brackets  with  subscript  e.  The  dissipative  part  is  the  remainder. 

Introduction  of  M is  purely  to  keep  the  expression  symmetrical. 

z 
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